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PACS 73 . 20 . -r - Electron states at surfaces and interfaces 

PACS 73 . 20 . At - Surface states, band structure, electron density of states 

PACS 73 . 21 . Ac - Multilayers 

PACS 81 . 05 .Uw - Carbon, diamond, graphite 

Abstract. - We report the existence of zero energy surface states localized at zigzag edges of Al- 
layer graphene. Working within the tight-binding approximation, and using the simplest nearest- 
neighbor model, we derive the analytic solution for the wavefunctions of these peculiar surface 
states. It is shown that zero energy edge states in multilayer graphene can be divided into three 
families: (i) states living only on a single plane, equivalent to surface states in monolayer graphene; 
(ii) states with finite amplitude over the two last, or the two first layers of the stack, equivalent to 
surface states in bilayer graphene; (iii) states with finite amplitude over three consecutive layers. 
Multilayer graphene edge states are shown to be robust to the inclusion of the next nearest- 
neighbor interlayer hopping. We generalize the edge state solution to the case of graphite steps 
with zigzag edges, and show that edge states measured through scanning tunneling microscopy 
and spectroscopy of graphite steps belong to family (i) or (ii) mentioned above, depending on the 
way the top layer is cut. 



Introduction. — In the past few years carbon physics 
' presented new challenges to the scientific community, in- 
. creasing the list of rather unusual phenomena occurring in 
' this life support element. On one hand, the discovery of 
. metal free carbon-based magnetism open a new research 
' field in fundamental physics, with possible applications 
in spin electronics [1-3]. On the other, the isolation of 
a single graphite layer - graphene - revealed an ultra- 
: relativistic system full of unconventional electronic prop- 
' erties, and regarded with great expectation from the point 
of view of applications [4-6] . 

The origin of the observed magnetism in carbon- 
based materials is still under debate, but the pres- 
ence of open edges seem to be an ubiquitous feature 
[3]. In proton bombarded graphite, which shows room 
temperature ferromagnetism, proton irradiation induces 
hydrogen-terminated edges [7, 8]. In activated carbon 
fibers and graphitized nanodiamond particles - known 
as nanographite - Curie- Weiss behavior and an enhanced 
paramagnetic susceptibility has been reported [2] . In these 
nanographites edges play a predominant role due to the 
built-in nano-dimension. Edges are assumed to induce tt- 



localized electrons due to surface (edge) states, which has 
been seen as a key ingredient to understand carbon's mag- 
netic behavior [1,3]. Indeed, the existence of edge states 
localized at zigzag edges of single layer graphene, induced 
either by extended defects or vacancies, is now well docu- 
mented and their magnetic behavior has been extensively 
reported [3,9-12]. 

Despite the positive correlation between edge state mag- 
netism in graphene single layer and magnetic phenomena 
in graphite and nanographite, strictly speaking, neither of 
them are a single layer of graphene. Although the inter- 
layer coupling is known to be very small, its effect is not 
negligible. To give an example, massless Dirac fermions 
in single layer graphene turn out to be massive in bilayer 
graphene [6] . This brings about the question whether edge 
states are robust to stacking, or in other words, whether 
multilayer graphene can support edge states localized on 
zigzag edges. Moreover, with the advent of graphene 
physics, also graphene multilayers (bilayer, trilayer, ...) 
were isolated. These graphene multilayers show interest- 
ing properties on their own [4,6], dissimilar from their 
single layer constituent, and can be even more suitable for 
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Fig. 1: (a) Side and top views of bilayer grapiiene, and its two 
families of edge states: monolayer (M) and bilayer (B); the 
vertical axes represent the associated charge densities (squared 
amplitudes), (b) Trilayer family of edge states (T) occurring 
in trilayer graphene (vertical axes represent charge densities), 
and all other possible edge states (in schematic view) . (c) Edge 
states in A'^-layer graphene (in schematic view). 



some device applications [13-15]. Therefore, the question 
whether multilayer graphene possesses edge state physics 
is of paramount importance. 

In this Letter we show that zero energy states localized 
at zigzag edges do exist in multilayer graphene. Using the 
simplest first nearest-neighbor tight-binding model, we de- 
rive the analytical expression for multilayer graphene edge 
states and show that their number is always equal to the 
number of layers occurring at the edge. The effect of sec- 
ond nearest-neighbor interlayer hopping is considered, and 
the robustness of multilayer graphene edge states is shown. 
Finally, we generalize the edge state solution to graphite 
steps, where experimental evidence for edge states has 
been widely reported [3] . The theoretical solution given in 
this Letter agrees well with experimental findings. Also, 
we predict that edge states in graphite steps should be 
seen in scanning tunneling microscopy (STM) even when 
the step occurs underneath the first graphite layer. 



ModeL — We model yli?-stacked multilayer graphene 
as shown in fig. 1(a) (for the simplest case of a bilayer), 
where non-interacting 7r-electrons are allowed to hop only 
between A and B sublattices. In what follows we use 
the terminology balcony layers for layers represented with 
dashed (red) lines, and non-balcony layers for those rep- 
resented with full (black) lines. Without loss of gener- 
ality we assume all edge atoms belong to the A sublat- 
tice. The zigzag edge breaks translational invariance along 
its perpendicular direction, enabling us to write an effec- 
tive one-dimensional Hamiltonian for a given momentum 
k € [0,27r[ along the edge (in units of a~^). The first 
nearest-neighbor tight-binding Hamiltonian can be writ- 



ten as 



i n 

-t± J2 + , (1) 



where ai-k,n ibi-k,n) is the annihilation operator at mo- 
mentum k and position n in sublattice Ai {Bi), i is the 
layer index and Dk = — 2cos(fc/2). The first term in 
eq. (1) describes in-plane hopping while the second term 
parametrizes the inter-layer coupling (t± <C t). The sym- 
bol • indicates a sum over non-balcony layers. Afterwards 
we consider the second nearest-neighbor interlayer hop- 
ping between A and B sublattices, which implies an extra 



term in eq. (1) given by -73 E„ 4fe,„(« 



-ik 



-ik/2jj 



kaizpi-k,n+i) + h-c, where 73 < t. 



Edge states in A^-layer graphene. — Multilayer 
graphene edge states are investigated by solving the 
Schrodinger equation, Hk IV'fe) = Ek \ipk)- The wavefunc- 
tion li/jk) is written as a linear combination of the site am- 
plitudes along the edge's perpendicular direction, \^k) = 
J2n \ai,k,n)+(3i{k,n) \bt,k,n)], where we have 

introduced the one-particle states \ci,k,n) — c|.^„|0), 
with Ci = ai,bi. In addition we require the boundary 
conditions ai{k,n 00) = ai{k^—l) = Pi{k,n — > 00) = 
Pi{k,—1) = 0, accounting for the existence of the edge at 
n = 0. Within our model, the Fermi energy of multilayer 
graphene always occurs at zero energy. Therefore, we ex- 
pect zero energy edge states to have interesting physical 
consequences, and we set Ek = 0- As a result, the two 
sublattices become completely decoupled, and only the 
sublattice to which edge atoms belong can support edge 
states.^ This means that we always have (3i{k,n) = 0. 

It was recently shown [17] that bilayer graphene sup- 
ports two types of zero energy edge states locahzed at 
zigzag edges for 27r/3 < k < An/ 3: one type restricted to 
the balcony layer and coined monolayer family, with am- 
plitudes equivalent to edge states in single layer graphene, 



a2{k,n) = a2(A;, 0)15^6" 



(2) 



and a new type coined bilayer family, with finite ampli- 
tudes over the two layers. 



ai{k, n) 
a2{k, n) 



ai(A:,0)D^e-*2" 



-ai{k,Q)Dl 



Dl 



(3) 



where the normalization constants are given by 
\a,{kM' = l-Dl and |ai(fc,0)|2 = (1 - Dir/[{1 - 



^In the ribbon geometry the two sublattices are equivalent, sup- 
porting edge states localized in opposite ribbon edges. In the semi- 
infinite system, only those localized in the edge sublattice survive 
[16,17]. 
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DI)"^ + t\/t'^]. The charge densities (squared amplitudes) 
associated with the two families of edge states arc repre- 
sented in fig. 1(a). Let us now consider a trilayer as shown 
in fig. 1(b), where a non-balcony layer is sandwiched 
between two balcony layers. Clearly, the bilayer family 
is not an edge state solution for this trilayer, as any 
finite amplitude at a non-balcony layer implies, through 
eq. (1), a finite amplitude over adjacent layers. We note, 
however, that our model ignores the coupling between 
next nearest-layers.^ Thus, if we construct a trilayer 
wavcfunction whose amplitudes over balcony /non-balcony 
layers mimic those for the bilayer family, it is guaranteed, 
apart from a normalization factor, that we have an edge 
state solution. More precisely, we arrive at a new type of 
edge state with finite amplitudes over three consecutive 
layers - trilayer family - whose analytic form can be 
written as 



tween two balcony ones, eq. (4). This is schematically rep- 
resented in fig. 1(c). Note that the munbcr of edge state 
families is always equal to the number of edge layers. 

Effect of 73. — In multilayer graphene the effect of 73 
is of the order of t±, and should be included in a consistent 
edge state solution. The buffer layer concept introduced 
previously, however, does not survive at a finite 73. In or- 
der to generalize the edge state solution to the present case 
we use the transfer matrix technique, following ref. [17]. 
For bilayer graphene the transfer matrix, defined as 
[aiik,n),a2ik,n)f = e-^W2T(2)"[ai(fc,0),a2(fc,0)]^, 
is given by 

'" u 



V 

Dk 



(5) 



T(2) = 

where u = Dk{l — v = 

X = -ifeJ''/^, with C = t^-fs/t^ 

pletely determined by the eigenvalues X± and eigenvectors 



_^g-ife/2(i _ £)2)^ and 

The edge states arc com- 



ai{k, n) 

0L2{k,n) 
a3{k,n) 



= —a2{k, O)!)^^^— e~'^^"~^^ (n — ) , of the transfer matrix. If |A±| < 1, then edge states ex- 



= a2(A;,0)£)^e-'^", 



-a2(fc,0)£) 



-s2)(n- 



(4) 



where the normalization constant is given by \a2{k, 0)^ = 
{l-Dlf /[{\-Dlf+2t\lt^]. The charge density (squared 
amplitude) associated with the trilayer family of edge 
states is represented in fig. 1(b). 

Additionally, the trilayer we have been discussing also 
supports edge states of the monolayer family localized at 
balcony layers, as schematically shown in fig. 1(b). In 
fact, this is a general result. Balcony layers have an edge 
sublattice which is not connected through t±_ to adjacent 
layers. Thus, the monolayer family is always an edge state 
solution in A'"-layer graphene. Even more generally, we can 
look at a balcony layer as a buffer layer. As can be seen 
from eq. (1), a finite amplitude over a balcony layer does 
not imply finite amplitudes over adjacent layers. For the 
trilayer shown at the bottom of fig. 1(b), where a bal- 
cony layer is sandwiched between two non-balcony layers, 
monolayer edge states certainly exist at the; micklle bal- 
cony layer. But because of the buffer layer character, also 
bilayer edge states are present, localized either at the two 
top or the two bottom layers. An immediate consequence 
of the buffer layer concept is the fact that the trilayer fam- 
ily is the most general edge state family we can have, and 
exists localized at any non-balcony layer and its two ad- 
jacent layers, with all other site amplitudes equal to zero. 
Therefore, we have three families of edge states occurring 
in multilayer graphene: (i) monolayer family for each bal- 
cony layer, eq. (2); (ii) bilayer family for each non-balcony 
layer that starts and/or ends the multilayer, eq. (3); (iii) 
trilayer family for each non-balcony layer sandwiched be- 



■^This is a reasonable approximation since in the Slonczewski- 
Weiss-McClure parametrization 72 , 75 <^tj_, 73 [6] . 



ist and are given by [ai{k,n),a2{k,v)Y' oc e 
apart from a normalization constant. Diagonalizing eq. (5) 
we obtain A± = £>fc(l - C/2) ± y/l^ D'H^/A - I) + 1. 
Simple algebra shows that for A+ the convergence condi- 
tion implies 2cos-^(VT+T/2) < A; < 2cos-i[-(l - C)/2] 
or 47r/3 < fc < 2 cos-^(-\/rTf/2), while for A_ it 
implies 2cos-'^{^/TTl/2) < k < 27r/3 or 2cos-^[(l - 
^)/2] < k < 2cos-i(-vTn/2). We conclude that bi- 
layer graphene still has two families of edge states for 
73 7^ 0, though the k range is slightly changed when 
compared with the 73 = case. In particular, we have 
only one family for k G [27r/3,2cos-^[(l - 0/2]] and 
k e [2cos-i[-(l - 0/2],47r/3], although the existence 
of edge states for k < 27r/3 and k > 47r/3 compensates 
this reduction, and we still have edge states for 1/3 of 
the possible A-'s, as in the 73 = 0. As a test to what 
has just been said, we have numerically computed the 
energy spectrum for a bilayer ribbon with zigzag edges 
{tj_ = 73 = 0.2t and width 400 tmit cells). The rcsuh 
is shown in fig. 2(a). Four flat bands at zero energy are 
clearly seen, and can be identified with the abovemen- 
tioned two families of edge states, two per edge. The insets 
reveal the k restrictions mentioned before. In fact, the val- 
ues of k that limit the existence or number of edge states 
coincide with the Dirac points and satellite Fermi points 
that arise when 73 ^ [18], as indicated by the thin red 
lines (the mismatch is due to the finite width of the ribbon, 
and consequent edge state overlap). The transfer matrix 
eigenvectors can be written as = [A± — D^, — ^e*'^/^]'", 
from which we can write two families of wavefunctions, 
\'ip±) = C± X]^o ^^^'^"^^^iX^j where the normalization 
constant is given by C± = [(1 — |A±|^)/(|A± — D^l^ -|- 
t\/t'^)]^/'^. Note, however, that and x~ are not 
orthogonal, implying the non-orthogonality of the two 
solutions \ip±)- It is convenient to orthogonalize \^-) 
with respect to jV'-i-)) whose result can be written as 
IV^-) = (IV-) - (V+|V^-)|V+))/(l - |(V'+IV-)n, where 
= C+C_{tl/t + iDl - 0/(1 -Dl+ 0- In 
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Fig. 2: (a) Energy spectrum for a bilayer ribbon with zigzag 
edges at finite 73. (b) Charge density for bilayer graphene edge 
states, monolayer (M) and bilayer (B) families, at fc/27r = 
0.35. Thin lines show the 73 = result. 



fig. 2(b) we show the squared amplitudes associated with 
\ip+) (left) and (right). The thin lines represent the 
edge states in bilayer graphene for 73 = 0, as given by 
eqs. (2) and (3). Clearly, as long as 73 <C we can 
identify with the monolayer family and {ip-) with 

the bilayer family discussed previously. For 73 ~ the 
edge state former monolayer family, already has an 

appreciable weight on both layers. 

The analysis made for bilayer graphene with 73 ^ 
can be extended to A^-layer graphene. Defin- 
ing the transfer matrix as [ai (fc, rt), . . . , Q;Ar(/c, n)]-'" — 
g-»fe»/2x(iV)»[ai(fc,0),...,ajv(A:,0)]^, we can derive its 
general pattern for a given number of layers N, 
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Dk 
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-^Dk 
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u 



T{N) 



(6) 

The periodic structure of T{N) is readily identified, and 
the transfer matrix for any multilayer graphene is eas- 
ily constructed. By diagonalizing eq. (6), and checking 
whether the eigenvalues A satisfy |A| < 1, we can conclude 
about the existence of edge states in A^-layer graphene for 
73 ^ 0. In fig. 3 the number of edge states per layer is 
shown in the plane number of layers (N) vs k. Panel 3(a) 
confirms what has been said for 73 = 0: same number of 
edge states as the number of layers for 27r/3 < k < 47r/3. 
As shown in panels 3(b) and 3(c), around the Dirac points 
the number of edge states for 73 7^ may be smaller than 
the number of layers, as we have seen for N = 2. How- 
ever, there is always a broad region in between the Dirac 
points where the number of edge states equals the num- 
ber of layers. Thus, we conclude that multilayer graphene 
edges states are robust to second nearest-neighbor inter- 
layer hopping. This result agrees with first principles cal- 
culations, where the presence of edge states and edge mag- 
netism have been reported in bilayer [19,20] and graphite 
[21] zigzag nanoribbons. 
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Fig. 3: Number of edge states per layer in A'^-layer graphene as 
a function of k: (a) 73 — 0; (b) 73 = t±; (c) 73 — 2t±. We set 
t± = 0.2t. 



Graphite steps. — Finally, we generalize the edge 
state solution to graphite steps, where experimental evi- 
dence for edge states has been widely reported [3,22-28]. 
The local density of states (LDOS) peak seen in STM 
of graphite steps has been interpreted as the experimen- 
tal confirmation of the theoretically predicted single layer 
edge states [29]. However, we can easily convince our- 
selves that monolayer edge states do not always provide 
an eigenstate for a zigzag step-edge. To see why, we con- 
sider fig. 4(a), where the two possible zigzag step-edges on 
the surface of graphite are shown. These two terminations 
are denoted a-type and /3-type. For an a-type termina- 
tion the edge carbon atoms occur exactly on top of carbon 
atoms of the underlying layer, while for a /3-type termi- 
nation they occur at the center of the hexagons. Obvi- 
ously, monolayer edge states as given by eq. (2) cannot be 
eigenstates for a-steps, as some finite amplitude must be 
induced on the second layer through eq. (1). For /?-steps, 
however, single layer edge states are indeed eigenstates. 

In order to have a step-edge we need at least two 
graphene layers. Indeed, step-edges are easily obtained 
from bilayer graphene just by growing one of the layers be- 
yond the edge. So, localized states at graphite steps can be 
understood by studying generalized bilayers where bottom 
and top graphene layers have different widths. The two 
families of edge states we have found to exist at zigzag 
edges of bilayer graphene remain eigenstates even when 
one of the layers is wider than the other. In particular, 
the solution given by eq. (2) and the non-orthogonalized 
solution for edge states of the bilayer family given by 



ai(fc, n) 
a2(fc,n) 



a,{k,0)D^e-^^\ 
-ai(fc,0)n-Driye~'*("-i^ 



(7) 



[a simple linear combination of eqs. (2) and (3)], can be 
adapted to the generalized bilayer just by adjusting the 
unit cell index n. Then, Gram-Schmidt orthogonalization 
gives the final solution. We should note, however, that 
the overlap between the two types of edge states is expo- 
nentially suppressed as the difference in layer width gets 
larger. When the width of one layer becomes infinite - the 
case of a perfect step - only one of the two possible solu- 
tions exists. Therefore, the possible localized solutions for 



p-4 



Surface states in multilayer graphene and graphite 




LDOS(P) - LDOS(a) 



G-©Y3 = 
□-□73 = 0.05? 



J \ I \ I \ L 



10 20 30 40 10 20 30 40 1 2 3 4 5 6 7 8 9 

n n n 




Fig. 4: (a) Possible zigzag steps on graphite's surface, 
(b) Charge density for edge states at a-type and /3-type step- 
edges {ka/2n — 0.35). (c) LDOS difference between /3- and 
Q-type steps as a function of n. (d)-(e) Top layer LDOS map 
for a- and /3-steps, respectively, (f)-(g) Underlying layer LDOS 
map for a- and /3-steps, respectively. We set t± = O.lt. 



zigzag step-edges (occurring at n — 0) are: 

atop(fc,n) = Cfei?^-^t", 
abottom(fc,") = -CA;7iL)^"-^-ie~' 

for an a- type step, and 



K«-i) 



(8) 



(9) 



for a /3-type step, where n > and the normalization con- 
stant in eq. (8) is given by \Ck\^ = (1 - Dl)'^/[{1 - Dlf + 
(1 -I- D\)t\/t^]. As in edge states discussed previously, 
the amplitudes in eqs. (8) and (9) refer to sites belonging 
to the same sublattice as the edge carbon atoms, while 
the amplitudes at the other sublattice are zero. Exam- 
ple charge densities for the given edge states are shown in 
fig. 4(b) for both the a-type and the /3-type step-edges. 

As can be seen from eqs. (8) and (9), or by inspection 
of fig. 4(b), there is an apparent asymmetry between the 
two families of edge states: edge states at /3-type steps 
live only on the top layer, while at a-type steps both the 
top layer and the underlying layer have a finite edge state 



amplitude. Consequently, we expect a similar asymmetry 
to be present in the LDOS peak induced by edge states 
at the Fermi level, which, ultimately, should be seen with 
STM. To better appreciate this effect, we have computed 
the LDOS of a generalized bilayer - bottom layer wider 
than the top layer - using the recursive Green's function 
method [30]. The calculated LDOS, which was accumu- 
lated in the range O.Olt near the Fermi energy, should be 
proportional, in the simplest approximation, to the local 
tunnel currents in the experimental STM images [25,31]. 
In fig. 4(c) we show, for the top layer (edge sublattice), the 
LDOS difference between /3-type and a-type terminations 
as we move away from the step at n = 0. As expected, 
the LDOS at /3-steps is higher and extends further into 
the bulk, a trend that is still present for realistic values 
of 73, as shown in fig. 4(c). This behavior agrees with 
STM results, where two types of edge states with different 
penetration depths have been seen [26]. Edge states with 
reduced penetration depth have been observed at a-type 
steps, whereas at /3-type steps the edge states extend fur- 
ther into the bulk, as we have obtained here for the top 
layer component. However, right at the edge, the STM 
intensity has been found to be higher at a-steps than [3- 
steps [26] . According to our analytical result the opposite 
should be seen. This discrepancy is most probably due to 
edge state admixture, as experimentally both a-type and 
/3-type steps coexist on the same step-edge. 

In fig. 4(d) and 4(e) we show the top layer LDOS map 
for a- and /3-steps, respectively. The former presents not 
only a reduced penetration depth, as previously discussed, 
but also higher intensity at sites connected to the underly- 
ing layer through t_L, as opposed to standard LDOS maps 
on the surface of bilayer graphene and graphite. This be- 
havior is characteristic of edge states at a-type steps, as 
given by eq. (8). fig. 4(f) and 4(g) show the underlying 
layer LDOS map for a- and /3-steps, respectively. As edge 
states at a-steps [eq. (8)] have a finite amplitude over the 
underlying layer, the LDOS map for this layer shows an in- 
creased intensity at and near the step [fig. 4(f)], although 
the lattice discontinuity only exists at the other layer. As a 
consequence, we expect a-steps to be detected in STM ex- 
periments even when they occur underneath the top layer. 
This feature is not seen in /3-type steps [fig. 4(g)]. 

Conclusions. — We have demonstrated the existence 
of zero energy states localized at zigzag edges of multilayer 
graphene and graphite steps. Stability to the presence of 
interlayer hopping 73 has been shown. The electron-hole 
symmetry breaking terms 74 (interlayer) and t' (inplane) 
are expected to induce edge state dispersion, but not to 
qualitatively modify the present results [32,33]. It should 
be noted that only perfect zigzag edges have been dis- 
cussed here. However, we expect edge state properties 
to be present in multilayer graphene and graphite steps 
even for irregular edges, as long as some zigzag units are 
present, as recently demonstrated for single layer graphene 
[34,35]. On the other hand, zigzag edges have been re- 
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cently observed in epitaxial graphene monolayer [36] , pro- 
viding the first indication tliat edge sliape can be a con- 
trollable parameter in the future. Our findings are rel- 
evant in the context of carbon based magnetism, where 
edge states seem to play an important role [1,3], and also 
in the context of graphene physics, where the reported 
self-doping in monolayer graphene [36] and suppression of 
conductance fluctuations near the neutrality point in bi- 
layer and trilayer graphene [37] can be seen as edge states 
driven effects. 

* * * 

E.V.C., N.M.R.P., and J.M.B.L.S. acknowledge 
financial support from POCI 2010 via project 
PTDC /FIS /64404 /2006. 

REFERENCES 

[1] Makarova T. and Palacio F., (Editors) Carbon Based 

Magnetism (Elsevier, Amsterdam) 2006. 
[2] Enoki T. and Kobayashi Y., J. Mater. Chem. , 15 

(2005) 3999 . 

[3] Enoki T., Kobayashi Y. and Fukui K., Int. Rev. Phys. 

Chem. , 26 (2007) 609. 
[4] Geim a. K. and Novoselov K. S., Nat. Mater. , 6 

(2007) 183. 

[5] Katsnelson M. I., Mater. Today , 10 (2007) 20. 

[6] Castro Neto A. H., Guinea F., Peres N. M. R., 

Novoselov K. S. and Geim A. K., The electronic prop- 
erties of graphene arXiv:0709.1163 (to appear in Rev. 
Mod. Phys.). 

[7] Kopelevich Y. and Esquinazi P., J. Low Temp. Phys. 

, 146 (2007) 629. 
[8] Ohldag H., Tyliszczak T., Hohne R., Spemann D., 

Esquinazi P., Ungureanu M. and Butz T., Phys. Rev. 

Lett. , 98 (2007) 187204. 
[9] Pereira V. M., Guinea F., Lopes dos Santos J. 

M. B., Peres N. M. R. and Castro Neto A. H., Phys. 

Rev. Lett. , 96 (2006) 036801. 
[10] Wakabayashi K., Electronic and magnetic properties 

of nanographite in Carbon Based Magnetism, edited by 

Makarova T. and Palacio F., (Elsevier, Amsterdam) 

2006 Ch. 12 pp. 279-304. 
[11] Lehtinen p. O., Foster A. S., Ma Y., Krashenin- 

NiKOV A. V. and Nieminen R. M., Phys. Rev. Lett. , 93 

(2004) 187202. 

[12] Son Y.-W., Cohen M. L. and Louie S. G., Nature , 

444 (2006) 347. 
[13] Castro E. V., Novoselov K. S., Morozov S. V., 

Peres N. M. R., Lopes dos Santos J. M. B., Nilsson 

J., Guinea F., Geim A. K. and Castro Neto A. H., 

Phys. Rev. Lett. , 99 (2007) 216802. 
[14] OosTiNGA J. B., Heersche H. B., Liu X., Morpurgo 

A. F. and Vandersypen L. M. K., Nat. Mater. , 7 (2008) 

151 . 

[15] Lin Y.-M. and AvouRis P., Nano Lett. , 8 (2008) 2119. 
[16] Wakabayashi K., Fujita M., Ajiki H. and Sigrist M., 
Phys. Rev. B , 59 (1999) 8271. 



[17] Castro E. V., Peres N. M. R., Lopes dos Santos J. 

M. B., Castro Neto A. H. and Guinea F., Phys. Rev. 

Lett. , 100 (2008) 026802. 
[18] McCann E. and Fal'ko V. I., Phys. Rev. Lett. , 96 

(2006) 086805. 

[19] Lee H., Son Y.-W., Park N., Han S. and Yu J., Phys. 

Rev. B , 72 (2005) 174431. 
[20] Sahu B., Min H., MacDonald A. H. and Banerjee 

S. K., Phys. Rev. B , 78 (2008) 045404. 
[21] Miyamoto Y., Nakada K. and Fujita M., Phys. Rev. 

B , 59 (1999) 9858 . 
[22] Klusek Z., Waqar Z., Denisov E. A., Kompaniets 

T. N., Makarenko I. v., TiTKOV A. N. and Bhatti 

A. S., Appl. Surf Sci. , 508 (2000) 161. 
[23] NiiMi Y., Matsui T., Kambara H., Tagami K., 

TSUKADA M. and Fukuyama H., Appl. Surf Set. , 241 

(2005) 43. 

[24] Kobayashi Y., ichi Fukui K., Enoki T., Kusakabe K. 
and Kaburagi Y., Phys. Rev. B , 71 (2005) 193406. 

[25] NiiMi Y., Matsui T., Kambara H., Tagami K., 
Tsukada M. and Fukuyama H., Phys. Rev. B , 73 

(2006) 085421. 

[26] Kobayashi Y., Fukui K., Enoki T. and Kusakabe K., 
Phys. Rev. B , 73 (2006) 125415. 

[27] Banerjee S., Sardar M., Gayatiiri N., Tyagi A. K. 
and Raj B., Appl. Phys. Lett. , 88 (2006) 062111. 

[28] SuGAWARA K., Sato T., Souma S., Takahashi T. and 
SUEMATSU H., Phys. Rev. B , 73 (2006) 045124. 

[29] Fujita M., Wakabayashi K., Nakada K. and Kusak- 
abe K., J. Phys. Soc. Jpn. , 65 (1996) 1920. 

[30] Haydock R., The recursive solution of the schrodinger 
equation in Solid State Physics, edited by Ehrenreich 
H., Seitz F. and Turnbull D., Vol. 35 (Academic Press, 
New York) 1980 p. 215. 

[31] Tersoff J. and Hamann D. R., Phys. Rev. B , 31 (1985) 
805 . 

[32] Peres N. M. R., Guinea F. and Castro Neto A. H., 

Phys. Rev. B , 73 (2006) 125411. 
[33] Sasaki K., Murakami S. and Saito R., Appl. Phys. 

Lett. , 88 (2006) 113110. 
[34] KuMAZAKi H. and Hirashima D. S., J. Phys. Soc. Jpn. 

, 77 (2008) 044705. 
[35] BiiowMiCK S. and Shenoy V. B., J. Chem. Phys. , 128 

(2008) 244717. 

[36] de Parga a. L. v., Calleja F., Borca B., Jr M. C. 

G. P., HiNAREjo J. J., Guinea F. and Miranda R., 

Phys. Rev. Lett. , 100 (2008) 056807. 
[37] Staley N. E., Puls C. and Liu Y., Phys. Rev. B , 77 

(2008) 155429. 



p-6 



